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Abstract 

Recent observational evidences of accelerating phase of the universe strongly 
demand that the dominating matter in the universe is in the form of dark en- 
ergy. In this work, we study the evolution of the apparent and event horizons 
<**] . for various dark energy models and examine their behavior across phantom 

^ ' barrier line. 
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1 INTRODUCTION 



The prediction of standard cosmology to have at present a phase of deceleration was 
ruled out in recent past by a series of observations namely the discovery of 16type la 
supernova(SNIa)by Riess et. al. (2004), WMAP(2003) and SDSS(2004). Using the Hubble 
telescope these observations has provided a distinct scenario of accelerated expansion of 
the present day universe. Thus a modification of Einstein equations [] becomes essential 
to incorporate this observational fact. One can either modify the geometry (i.e., the left 
^ ' hand side of Einstein equation ) or the matter itself(i.e., the R.H.S.) if not both. Due to 

£ | modification of geometry, one can introduce modified gravitytheory namely f(R) gravity, 

Brane scenario etc while change in the matter part indicates inclusion of some unknown 
kind of matters having large negative pressure so that strong energy condition (p+3p > 0) 
is violated. Such an unknown matter is known as dark energy(DE). 
■ In literature, there are various DE models to match with observational data. The 

simplest model representing DE is the Cosmological Constant which was introduced by 
Einstein himself, surprisingly many years before the starting of DE craze. However, this 
model of DE is not very popular due to many inherent drawbacks (for example fine 
tuning problem ( Steinhardt 1997)). The other candidates for DE are variable cosmological 
constant (Shapiro et. al. 2009; Sola et al. 2005 ; Solaet al. 2006 ), the canonical scalar 
^ ■ field (Dutta et. al. 2009;Guo et al. 2007; Liddle et. al. 1999;Ratra et al. 1988;Wetterich 

, 1988; Zlatev et al. 1999) (quintessence field), scalar field with negative kinetic energy 

(phantom field) (Caldwell 2002 ; Caldwell et al. 2003; Nojiri et al. 2003B ; Onemli et al. 
2004; Saridakis 2009; Setare et al. 2008; Setare et al. 2009) or a quintom field (Capozziello 
et al. 2006; Elizalde et al. 2004; Feng et al. 2005; Feng et al. 2006 ; Guo et al. 2005 ; Li 
et al. 2005 ; Setare 2006 ; Setare et al. 2008 A ; Setare et al. 2008B; Setare et al. 2008C 
; Setare et al. 2009A; Zhao et al. 2006 ;) (a unified model of quintessence and phantom 
field). Further a combined effort of quantum field theory and gravity leads to speculate 
some nature of DE and is known as holographic dark energy (HDE) model ( Copeland 
et al. 2006; Durrer et al. 2008; Nojiri et al., 2007; Padmanabhan 2002; Sahni 2005 ,2006; 
Nojiri et. al. 2006B). 

In the present work we study the evolution of the horizons(apparent and event) for 
different DE models namely (a) DE with barotropic equationof state, (b) holographic 
DE(HDE) and (c) a non interacting two fluid system-HDE and dark matter in the form 
of dust. The paper is assigned as follows : Basic equations are presented in the section 
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2, evolution of the horizons are studied for the above three matter systems in section 3, 
section 4 deals with thermodynamical analysis of the universe bounded by the horizons. 
The paper ends withdiscussion and concluding remarks in section 5. 



2 Basic equations 

For simplicity let us start with homogeneous and isotropic model of the universe (namely 
Friedmann- Robertson- Walker(FRW) model), having line element 



ds 2 = —dt 2 + a 2 (t) 



1 — kr 2 
= h aa dx a dx" + R 2 dQ, 



dr2 , 2 , n 2 

+ r dil 



(1) 



where 



and 



h a t = diag [ — 1, — — — - ) , (a, b = 0, 1 with x° = t, x 1 = r) 
V 1 — kr 2 J 

dfl 2 = d9 2 + sin 2 9d(j} 2 is the metric on unit two sphere. 



R = ar is the radius of the sphere(area-radius) , 'a' is the scale factor and k = 0, ±1 stands 
for flat, closed and open model of our universe respectively. 

The matter is chosen as a perfect fluid with energy momentum tensor 

T,j, v = (p+p) u^Uy — pg^ v (2) 

So the Einstein field equations are (choosing 8ttG = 1 = c) 



H-^ = -\{p+V) (4) 
a 2 2 

and the energy conservation equation is 

p + 3H(p + p) = (5) 

Combining Q and (j4| we get, 

^ + # 2 = ^ = -i(p + 3 P ) (6) 

The dynamical apparent horizon which is essentially the marginally trapped surface 
with vanishing expansion, is defined as a sphere of radius R = Ra such that 

h ab d a Rd b R = (7) 

which on simplification gives 

Ra = 1 (8) 
lH 2 + 4 



The event horizon on the other hand is defined as (Davis 1998) 
Re = —a sinh(r) k = —1 
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Re = -ar k = (9) 



Re = —a sin(r) k — +1 

where r is the usual conformal time defined as 



\t\ < oo (10) 



dt 

Note that if |r| = oo, event horizon does not exist. Also the Hubble horizon is given 

by 

Rb = jj (11) 
The horizons are related by the following relations (Mazumder 2009): 

Ra=Rh < Re forK = 

Rh<Ra<Re forK = -l 
Ra < Re < Rh 

or forK = +1 

Ra < Rh < Re 

3 Evolution of the horizons and consequences 

The time variation of the horizon radii are given by 

Ra =-h(h- ^) R\ (12) 



lh = II Hl - \/l-^R% (13) 



Rh = (14) 



~H~ 2 

One may note that the expression for Re given in references (Davis 1998) and 
(Mohseni Sadjadi 2006) are true only for k = 0. So the theorems given in the papers 
of Davis(1998) and Sadjadi(2006) are only valid flat universe. However, in the present 
work from the above expression (i.e., equation (13)) we see that Re is an increasing or 
decreasing function of time that depends only on whether Re > or < Ra- it does not 
depend on the nature of the matter involved as claimed by Davis and Sadjadi. 

We shall now study the variation of the horizons with the evaluation of the universe. 
Due to observed accelerating phase of the universe, the matter is assumed to be in the 
form of the DE having equation of state p = uip. 

Case I : uj is constant 

For simplicity, if we assume the flat model of the universe then from equation © we have 
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with a = | (1 + co). Hence in the quintessence era we have < a < 1. Now solving the 
Einstein field equation © and the conservation equation © we have 



a — aot", p — pot 2 (15) 



Then the horizons are given by 



nt 

Re=, , R A = at (16) 

1 — a 

Hence, over one Hubble time (tn = jj) both have the same time variation, i.e., 

tn^=a (h = EorA) (17) 
Thus there are no significant changes of the two horizons over the Hubble time. 

Case II : u) is variable 

Here the choice of DE is holographic model. The holographic principle states that the no. 
of degrees of freedom for a system within a finite region should be finite and is bounded 
roughly by the area of its boundary. From the effective quantum field theory one obtains 
the Holographic energy density as (Cohen et al. 1999) 

3c 2 
He 

where free dimensionless parameter c is estimated from observation and IR cut off is 
chosen as Re to get correct expression. 

Then using expression (|18l) in the conservation equation ((5} the expression for the 
equation of state parameter is given by 



— HVt-* < I9) 

where fio = and Qk = -^r^i are the density parameters corresponding to DE 

and curvature respectively. Now from equations (1 12[) and (|13p the time variation of the 
horizons over one Hubble time are given by 



te ^ = |(i + w) = i-J^-n fc (20) 



and 



••E- 1 -^-"" (21) 

So both the horizons have the same time variation over one Hubble time. 

Case III : Variable u and two fluid syatem 

Here we consider a non-interacting two fluid system having one component in the form 

of HDE and the other component as dark matter (in the form of dust). So the Einstein 
equations for flat FRW model now become 

H 2 = ^(p D +p m ) (22) 
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H = ~(po + Pm+PD) (23) 
As the fluids are non-interacting so the energy conservation equations are 

p m + 3H Pm = (24) 

and 

p D +3Hp D (l+u) =0 (25) 

So from the expression of the energy density for the HDE(given by equation Q18|ll we 
have as before 

i 2 VnE ( . 

" = ~3"3^ (26) 
where variation of the density parameter is given by 

«' = fif,(l-fi D ){-L + _J=) (27) 

Now the change of the horizons over one Hubble time are given by the expressions 

'■k = I " 5 fi ° - (28) 

Thus compared to between the the above two cases the changes of the horizons over one 
Hubble time are not identical, though they do not change significantly. 



4 Thermodynamics of the Universe and the role 
of the horizons : 

Here we consider the universe bounded by the event or apparent horizon as a thermody- 
namical system. In the previous section we have shown that neither the apparent horizon 
nor the event horizon change significantly over one Hubble time scale so equilibrium ther- 
modynamics can be applied here with temperature and entropy on the horizon similar to 
black holes. 

Case I : Matter in the form of perfect fluid : 

Here matter bounded by the horizon is considered to be in the perfect fluid. The total 
entropy change can be written as (for details see Mazumder et al 2009) 

^ {Sl + S h ) = ^(p + p)R h (30) 

where Rh is the radius of the horizon(event or apparent), Si and Sh are respectively the 
entropy of the matter bounded by the horizon and that of the horizon, p and p are the 
energy density and the thermodynamic pressure of the inside matter and Tf, is the temper- 
ature of the horizon as well as of the inside matter for equilibrium thermodynamics. Thus 
generalised second law of thermodynamics will be valid in quintessence era (p + p > 0) if 
the radius of the horizon increases with time while in phantom era (p + p < 0) the radius 
of the horizon should decrease. 
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Case II : Matter in the form of HDE : 

If we differentiate the expression for the energy density of the HDE (i.e., equation (JTSJl) 
then using the energy conservation equation |[5J we obtain (after simplification) 

Re = ^HR E (l+u) (31) 
As before variation of the total entropy is given by equation (|30[1 which using (|31[) becomes 

i(S I + S E )^ b ^ PD (l + uf (32) 

for event horizon. For the apparent horizon using equation (j 12|) and the Friedmann equa- 
tion @, equation (|30p simplifies to 

^ {Sl + SA)= 2 J^k pD{1 + u f (33 ) 

Thus generalised second law of thermodynamics hold for both the horizons when matter 
is purely in the form of HDE. 

Case III : Non-interacting two fluid system : 

Here matter in the universe bounded by the horizon (event or apparent) is in the form of 
non-interacting two fluid system-one component is HDE (pd, Pd) and the other is dark 
matter in the form of dust (p m ). Then total entropy variation (for details see Mazumder 
et. al. 2010) is given by 

£(S I + S A ) = ^^{p m + p D (l + u } )}R h (34) 
at 1 h 

Thus energy density of dark matter plays a key role for the validity of the generalized 
second law of thermodynamics particularly in phantom era. 



5 Discussions and Concluding remarks: 

We shall now discuss the behavior of the horizons with the evolution of the universe both 
in Quintessence and Phantom eras. From the conservation equation (5) we see that in 
Quintessence era p is monotonic decreasing which reaches a local minima at the phantom 
crossing and increases again with the evolution of the universe as shown in Fig I. So the 
matter density has some short of bouncing behavior at the phantom crossing. However, 
if the universe starts contracting in phantom era (i.e., H < 0) then conservation equation 
demands p should still decreases in the phantom era and there is a point of inflexion at 
the phantom barrier as shown in Fig. II . For both the possibilities in phantom era p has 
peculiar behavior when matter is exotic in nature (i.e., p + p < 0). In the first case when 
universe is expanding p also increases in the phantom era indicating some matter creation 
phenomena (of unknown nature) during that epoch. On the other hand, when universe 
starts contraction in the phantom era, p still decreases, indicating destruction of mass in 
that era. 

We shall now present a comparative study of the evolution of the horizons across the 
phantom barrier with the expansion of the universe for the various DE matter distribu- 
tion in the following tabular form 



G 



Energy Density Curve 
For 

Expanding Universe 




Fig.I 

Fig. I represents the variation of energy density with the evolution of the universe 
in an expanding model. The dotted vertical line denotes the phantom divide or 
phantom barrier line. 
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Fig.III 

Fig. Ill represents variation of event horizon and the apparent horizon respec- 
tively in an expanding universe model. The dotted vertical line again denotes the 
phantom divide line. As the previous diagrams left side of which is denoting the 
quintessence era whereas the right hand side represents the phantom era. 
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Quintessence Era 


Phantom Era 



1 



Fig.IV 

Like Fig.III, Fig.IV also represents two curves showing the variation of event hori- 
zon and apparent horizon. The dotted vertical line denotes the phantom divide 
line. The left side of which is denoting the quintessence era whereas the right hand 
side represents the phantom era. 
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Table I:comparative study of the two horizons across phantom barrier 





Quintessence Era 


Phantom Era 


Apparent 
Horizon 


(a) For HDE model apparent horizon 
increases with the evolution of the universe. 

(b) For other DE model the behavior 

is identical as HDE. 


(a) For HDE model Ra decreases with time. 

(b) Ra has same behavior for other 
DE model as in HDE. 


Event 
Horizon 


(a)For HDE model event horizon 
increases and hence from equation(|13p 
we always have Re > Ra- 
(b) Re has similar behavior as in HDE, (for 
other DE models). 


(a) For HDE model Re decreases with time 
and so we have Re < Ra- The variations 
are schematically shown in figure III 
(b) For other DE model, Re may have 
behavior as in HDE model. Also Re may 

still be an increasing function with the 
evolution and the nature is schematically 
shown in figure IV. 



We have also studied in the last section the thermodynamics of the universe bounded 
by the horizons with different DE models. We have assumed the validity of the first law 
of thermodynamics and examined whether the GSLT holds or not. The conclusions from 
this thermodynamical study has been presented below in table II. 

Table II: Validity of GSLT for different DE model 





Quintessence Era 


Phantom Era 


Universe 
bounded 

by the 
apparent 

horizon 


(a)For HDE model GSLT is always 
satisfied. 

(b) For other DE model GSLT is satisfied 
throughout the evolution, 
(c) For non-interacting 2-fluid system 
there is always validity of GSLT. 


(a) GSLT is always satisfied for the HDE model. 

(b) Across the phantom barrier GSLT is satisfied. 

(c) GSLT will be satisfied provided |1 + w| > 


Universe 
bounded 
by 
the 
event 
horizon 


(a) Thermodynamical system respects 
GSLT for HDE model. 
(b)GSLT is obeyed for other DE model, 
(c) No restriction is needed for the validity 
of GSLT for noninteracting 2-fluid system. 


(a) No restriction is needed for the validity 
of GSLT. 

(b) Validity of GSLT depends on the behavior 
of Re- GSLT will be respected for the 
fig III while it will be violated for fig. IV. 
(c) GSLT will be satisfied for variation 
of Re according to fig III provided 
|l+w| < 

1 Pn 



Thus, from the above study we see that both the evolution of the horizons as well as 
the matter density have some strange behavior in the phantom era, i.e., across the phan- 
tom barrier line, therefore, for future work cosmological evolution in phantom region will 
be done more in details and also it will be interesting to explain the particle creation 
in the phantom era with the expansion of the universe and possibly the mechanism of 
particle creation may remove the possible future singularity. 
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